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求解 0-1 背包 问题 的 混合 蝙蝠 算法 
万 晓 琼 ， 张 惠珍 、 
(上 海 理工 大 学 管理 学 院 ， 上 海 200093) 


摘 要 : 针对 基本 蝙蝠 算法 易 陷 入 局 部 最 优 、 收 你 速度 慢 等 缺点 ， 对 其 进行 优化 研究 。 基 于 0-1 背包 问题 的 具体 特征 ， 
在 基本 蝙蝠 算法 原 有 概念 和 框架 的 基础 上 ， 引 入 遗传 算法 中 的 交 又 机 制 以 及 反 置 算 子 建立 全 新 的 位 置 转移 方式 和 局 部 
搜索 规则 ; 加 入 贪心 策略 进行 解 的 可 行 化 和 充分 利用 ， 增 强 局 部 搜索 能 力 ， 加 快 算法 收 化 速度 ， 构 建 全 新 的 混合 蝙 晤 
算法 。 将 混合 蝙蝠 算法 应 用 于 两 组 0-1 背包 算 例 ， 仿 真实 验 结果 优 于 自 适应 元 胞 粒子 群 算 法 、 基 本 蝙蝠 算法 和 贪心 二 
进 制 蝙蝠 算法 。 结 果 验 证 了 该 混合 算法 求解 0-1 背包 问题 的 可 行 性 和 有 效 性 。 
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Hybrid bat algorithm for solving 0-1 knapsack problem 


Wan Xiaoqiong, Zhang Huizheni 
(School of Management, University of Shanghai for Science & Technology, Shanghai 200093, China) 


Abstract: In order to overcome the shortcomings of the basic bat algorithm, such as easy to fall into local optimum and slow 
convergence speed, the optimization study was carried out. This paper proposed a new hybrid bat algorithm based on the original 
features, framework of the basic bat algorithm and specific characteristics of the 0-1 knapsack problem. The hybrid algorithm 
introduced the crossover mechanism of the genetic algorithm and the inverse operator to construct a brand-new location transfer 
method and local search rule. The greedy strategy was added to make the solution feasible and fully utilized, enhancing the local 
search ability and speeding up the convergence. This hybrid bat algorithm was applied to two sets of 0-1 knapsack cases. The 
simulation results exceeded adaptive cellular particle swarm optimization, basic bat algorithm and greedy binary bat algorithm. 


The results verify the feasibility and effectiveness of the hybrid algorithm in solving 0-1 knapsack problem. 
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0 ”引言 为 两 类 : 精确 算法 和 启发 式 算法 。 精确 算法 如 : 分 支 定 界 法 ?9、 
动态 规划 法 %、 回 溯 法 %、 贪心 算法 "等 。 精确 算法 虽然 可 以 求 得 

背包 问题 (knapsack problem, KP) 是 运筹 学 中 一 个 典型 的 优化 问题 的 最 优 解 ， 但 无 法 解决 随 着 物品 数量 的 增加 而 产生 的 
NP-Complete 组 合 优 化 问题 ， 则 在 寻求 满足 背包 约束 的 条 件 下 “组 合 爆炸 ”问题 ， 启 发 式 算法 的 出 现 ， 为 解决 此 类 问题 开辟 
有 具 有 最 大 价值 的 物品 装载 方案 ， 现 实生 活 中 的 诸多 问题 都 可 被 了 一 条 新 途径 。 现 在 已 有 多 种 启发 式 算法 被 用 来 解决 0-1 背包 
归结 为 背包 算 例 ， 如 投资 决策 问题 、 任 务 调度 问题 、 货 物 装载 问题 ， 如 蚁 群 算法 ””、 和 声 算法 "、 遗 传 算 法 、 粒 子 群 算 
] 题 、 材 料 切 制 问题 等 。 背包 问题 一 般 可 以 分 为 三 类 : 0-1 背包 法 “" 、 狼 群 算法 路、 鲸鱼 算法 "、 晴 蜂 算 法 "和 等。 上述 算 法 在 
题 (0-1 knapsack problem, 0-1 KP) 中 、 完 全 背包 问题 和 多 解决 0-1 背包 问题 时 取得 了 一 定 的 成 果 , 同时 也 显现 出 易 早 熟 、 
背包 问题 ， 其 中 0-1 背包 问题 是 最 基本 的 背包 问题 ,包含 了 背 后 期 收敛 速度 慢 等 不 足 。 为 了 克服 单一 启发 式 算法 的 缺点 ， 运 
题 设计 状态 以 及 方程 的 最 基本 思想 ， 其 他 类 型 的 背包 问题 用 混合 算法 求解 组 合 优化 问题 成 为 趋势 ， 如 混合 遗传 算法 、 
以 转化 为 0-1 背包 问题 进行 求解 。 本 文 则 在 对 0-1 背包 问 基于 遗传 算 子 的 蚁 群 算 法 ”"、 基 于 贪心 策略 的 遗传 算法 "等 。 
进行 研究 。 坊 蝠 算法 (bat algorithm, BA) “是 剑桥 大 学 学 者 Yang 
0-1 背包 问题 自 提 出 以 来 , 一 直 是 众多 学 者 研究 的 热点 。 从 ”Xinshe 基于 蝙蝠 回声 定位 行为 于 2010 年 提出 的 一 种 启发 式 算 
已 有 的 研究 成 果 可 以 看 出 ， 目 前 求解 该 问题 的 算法 主要 可 以 分 法 ， 具 有 参数 少 、 结 构 简 单 、 鲁 棒 性 强 等 优点 。 然 而 ， 与 其 他 
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启发 式 算 法 相似 ， 基 本 蝙蝠 算法 也 存在 易 陷 入 局 部 最 优 和 路 鲍 蝙蝠 算法 。 

速度 慢 等 缺点 。 针 对 这 些 缺 点 ， 本 文 在 基本 蝙 蝠 算法 原 有 概念 。 2.1.1 初始 化 蝙蝠 位 置 规则 

和 结构 的 基础 上 ， 引 入 遗传 算法 中 的 交叉 机 制 提 高 全 局 搜索 能 根据 0-1 背包 问题 解 的 性 质 ， 本 文采 用 二 进 制 编码 方式 随 
力 ， 加 入 反 置 算 子 和 贪心 策略 增强 局 部 搜索 能 力 ， 设 计 出 全 新 机 初始 化 蝙蝠 种 群 位 置 。 用 1 *n 维 向 量 xi = (cia， Xiz Xin) 
的 混合 蝙蝠 算法 (hybrid bat algorithm, HBA ) 来 求解 0-1 背包  ” 表示 第 i 只 蝙蝠 的 初始 位 置 , 其 中 xij E {0,1},1<j <n, 当 
问题 ， 并 与 文献 [23,24] 的 0-1 背包 算 例 相对 比 ， 结 果 验 证 了 该 xij = 1 时 ,物品 j 装 入 背包 ， 反 之 ， 物 品 j 不 装 入 背包 。 具 
算法 在 寻 优 能 力 和 收敛 速度 上 均 有 较 大 优势 。 体 映 射 关 系 如 式 (4) 所 示 。 


1, rand > 0.5 
1 ”0-1 背包 问题 | 
百 包 门 是 0， rand < 0.5 0% 


0-1 背包 问题 是 经 典 的 组 合 优化 问题 , 具体 可 以 描述 为 : 给 2.1.2 蝙蝠 速度 与 位 置 更 新 规则 
定 有 重量 限制 的 背包 以 及 具有 重量 和 价值 的 nn 个 待 装 物 品 ， 求 1) 速度 更 新 规则 
解 目 标 为 在 背包 载重 限制 下 选择 物品 装 入 背包 并 使 得 装 入 物品 由 于 离散 优化 问题 的 特殊 性 , 在 初始 化 蜗 蝠 位 置 以 及 后 续 
价值 最 大 。 该 问题 数学 模型 表示 如 下 : 速度 定义 、 位 置 转移 编码 都 是 一 步 到 位 的 ， 故 删 去 基本 蝙 晤 算 
C 法 中 频率 三 这 个 参数 。 重 新 定义 第 只 蝙 蝙 x; 和 最 优 蝙蝠 x 
max F(x)= Di Xi (1) es ee ee i 
全 之 间 的 Manhattan 距离 为 蝙蝠 的 速度 。 则 第 i 只 蝙蝠 在 定义 的 
= 搜索 空间 中 第 t 时 刻 vt 的 更 新 公式 为 
s.t. > wixi: <C (2) n 
- vf = > ly 一 xj (5) 
口 ,4 EI E 
| Lg 物品 i 装 入 背包 中 viel (3) 人 
0， 否则 2) 位 置 更 新 规则 


模型 中 ,1 = {i|i= 1,2，...n} 为 物品 集 ;v = {vi|i= 本 文采 用 遗传 算法 中 的 交叉 机 制 进行 位 置 转移 , 一 只 蝙蝠 
1,2，...n}) 为 物品 价值 集 ; w = {wi | i= 1,2，...n} 为 物品 对 应 一 条 染色 体 ， 位 置 更 新 过 程 主要 可 以 分 为 两 步 : 

量 集 ; C 为 背包 载重 量 ; 目标 函数 式 (1) 表示 最 大 化 装 入 背 a) 子 代 遗 传 父 代 的 相同 基因 。 具 体 为 : 设 最 优 蝙 蝠 x, 为 父 
包 物 品 价值 ;约束 式 (2) 保证 装 入 背包 物品 的 总 重量 不 超出 代 1， 第 i 只 蝙蝠 xj 为 父 代 2， 对 比 父 代 1、2 相应 的 基因 位 ， 
背包 载重 量 ， 约束 式 (3) 中 xi 为 决策 变量 。 当 基 因 相 同 〈( 即 相应 基因 位 上 同 为 1 或 0) 时 ， 将 相应 基因 直 

i 接 遗 传 给 子 代 ; 当 基 因 不 同时 , 子 代 相应 基因 位 以 通配符 * 填 
nn 充 ， 记 此 时 的 子 代 蝙 同位 置 为 xyia。 
蝙蝠 算法 是 自然 界 中 蝙蝠 群体 回声 定位 行为 的 智能 模拟 ， b) xmia 中 通配符 的 个 数 为 vf， 即 为 第 i 只 蝙蝠 的 转移 速度 ， 
是 一 种 搜索 全 局 近似 解 的 智能 优化 算法 。 蝙 蝠 个 体 是 蝙蝠 算法 ”对 于 这 vt 个 不 同 的 基因 , 设置 追随 概率 pyo1, 建立 转移 公式 如 式 
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的 基本 单元 ， 它 们 通过 嘴巴 和 鼻孔 向 外 发 射 很 强 的 超声 波 ， 并 。”“〔6) 所 示 。 
通过 超凡 的 大 耳 廓 接收 回音 ， 再 依据 发 出 超声 波 和 接收 到 回声 ny = td (uF, prar) (6) 
之 间 的 时 差 、 到 达 两 耳 之 间 的 强度 差 、 回 声 定位 波形 变化 建立 体 为 ， 当 rand > proi 时 ， 子 代 追 随 最 优 蝙 蝠 ， 遗 传 父 代 


空间 三 维 场景 来 感知 猎物 或 障碍 物 的 距离 、 方 位 以 及 性 质 ， 进 1 的 对 应 基因 ， 反 之 ， 遗 传 父 代 2 的 对 应 基因 ， 最 后 更 新 子 代 
而 制定 捕捉 猎物 和 躲避 障碍 物 的 计划 。 算 法 思想 可 以 概括 为 : 蚁 蝠 位 置 为 xnew。 利 用 这 种 交叉 机 制 ， 一 方面 子 代 可 以 遗传 父 
蝙蝠 种 群 随 机 散布 于 n 维 问题 空间 内 ， 把 它们 映射 为 优化 问题 。 代 优 秀 的 基因 。 另 一 方面 , 加 入 追随 概率 pyo1， 增添 跟从 最 优 蝙 
空间 内 的 s 个 解 ， 所 优化 问题 的 目标 函数 对 应 解 的 适应 度 值 ， 蝠 的 随机 性 以 维持 蝙蝠 种 群 的 多 样 化 ， 一 定 程度 上 克服 了 基本 
适应 度 值 决定 解 的 优 务 ， 算 法 的 优化 过 程 即 为 蝙蝠 种 群 的 其 他 ”蝙蝠 算法 容易 早熟 的 缺点 。 以 n = 6 为 例 , 具体 交叉 重组 过 程 如 
个 体 调整 响 度 、 脉 冲 发 射 率 在 解 空 间 中 追随 最 优 蝙 蝠 的 迭代 过 图 1 所 示 。 


程 。 2.1.3 局 部 搜索 规则 

2.1 求解 0-1 背包 问题 的 混合 蝙蝠 算法 本 文采 用 的 局 部 搜索 方式 有 两 种 ， 一 为 利用 反 置 算 子 对 最 
蝙蝠 算法 最 初 被 用 来 求解 函数 优化 问题 ， 近 几 年 有 学 者 提 优 蝙 晤 进行 局 部 搜索 ;二 为 使 用 贪心 策略 对 不 可 行 解 和 非 充分 

供 新 的 改进 思想 和 方法 去 解决 离散 问题 。 在 本 文中 ， 为 了 把 蝙 。 利用 解 进行 局 部 搜索 。 

算法 应 用 于 0-1 背包 问题 ， 在 求解 优化 问题 的 过 程 中 ， 利 用 1) 反 置 算 子 

二 进 制 编码 方式 初始 化 蝙蝠 位 置 ， 在 蝙蝠 位 置 更 新 规则 中 引入 相对 于 其 他 智能 优化 算法 ， 蝙 蝠 算法 可 以 通过 比较 脉冲 发 

遗传 算法 的 交叉 机 制 ,利用 反 置 算 子 构建 全 新 的 局 部 搜索 规则 ， ” 生 率 7; 和 随机 数 的 大 小 ， 实 现 对 全 局 搜索 和 局 部 搜索 之 间 动 态 


称 之 为 混合 蝙蝠 算法 ， 下 文 从 初始 化 蝙蝠 位 置 、 蝙 蚁 速度 与 位 。 转换 。 即 当 rand > 7i 时 , 在 当前 最 优 解 集中 选取 一 个 最 优 蝙蝠 ， 
置 更 新 、 局 部 搜索 以 及 音量 和 脉冲 的 更 新 四 个 方面 来 介绍 混合 。 进行 随机 游 走 ， 产 生 局 部 新 解 。 由 于 二 进 制 编码 的 特殊 性 ， 基 
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本 蝙蝠 算法 的 局 部 搜索 方式 并 不 能 取得 优异 结果 。 本 文 引入 反 
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而 趋 于 相似 的 蝙蝠 种 群 容易 跳出 局 部 最 优 。 对 于 选 定 的 基因 位 


置 算 子 设置 全 新 的 局 部 搜索 规则 ， 具 体 为 : 反 置 概率 pnor 随机 。 jnor € Jnot 具体 优化 方法 如 下 : 
选取 最 优 蝙 蝠 的 mn* pnot 个 基因 进行 反 置 处 理 ,最 后 更 新 局 部 新 1， xj = 本 
% 9 一 
解 xnew。 这 种 随机 反 置 的 局 部 搜索 方式 使 得 随 着 先 代 次 数 增加 0, xj =1 
父 代 1 : 1 0 1 1 0 1 
父 代 2 : 1 1 1 0 1 1 
第 一 步 : 子 代 遗传 父 代 相同 的 基 
子 代 xmia 1|*|1|*|*|1 
第 二 步 : 当 rand > prol 时 ， 追 随 最 优 蝙 晤 
子 代 xnew 1 0 1 0 0 1 
图 1 蝙蝠 算法 位 置 更 新 过 程 


2) 贪心 策略 

a) 对 不 可 行 解 进行 贪心 策略 。0-1 背包 问题 作为 具有 约束 
条 件 的 组 合 优化 问题 ， 在 蝙蝠 位 置 随机 初始 化 、 位 置 更 新 以 及 
反 置 算 子 局 部 搜索 部 分 均 有 可 能 产生 无 效 解 ， 即 装 入 物品 总 重 
量 超出 背包 重量 限制 ， 针 对 这 些 不 可 行 解 若 直接 舍弃 ， 则 会 丢 
失 已 搜索 到 的 信息 ， 导 致 算法 收敛 速度 变 慢 ， 故 本 文 引入 贪心 
策略 对 不 可 行 解 进行 局 部 搜索 : 计算 已 装 入 背包 物品 的 价值 比 
依次 铭 去 低 价值 比 的 物品 ， 直 至 剩余 物品 总 重 
不 超出 背包 载重 。 

b) 对 可 行 解 进行 贪心 策略 。 由 于 初始 化 解 的 随机 性 和 可 行 
化 解 的 限制 ， 可 能 导致 背包 容量 的 不 充分 利用 ， 针 对 这 部 分 解 
采用 贪心 策略 进行 局 部 搜索 : 计算 未 装 入 背包 的 物品 价值 比 
voli = Vi/wi， 依 次 加 入 高 价值 比 的 物品 ， 直 至 背包 利用 率 达 到 
最 大 。 
2.1.4 蝙 蝙 音量 和 脉冲 发 生 率 更 新 规则 
在 蝙蝠 接近 猎物 时 ， 超 声波 音量 逐渐 降低 ， 脉 冲 发 生 率 逐 


加 


voli = Vi/wi， 


步 提 高 。 音 量 入 和 脉冲 发 生 率 x 的 更 新 公式 为 
At = aAt! (8) 
rt = {1 exp[—y(t — 1)]} (2) 


其 中 : a 为 音量 衰减 系数 , y 为 脉冲 发 生 率 增长 系数 ,，r? 为 最 大 

脉冲 频率 。 对 任意 0 < a < 1,y >0, 可 以 看 出 : 

rn, as t—»%。 

2.2 求解 0-1 背包 问题 的 混合 蝙蝠 算法 设计 
通过 对 基本 蝙蝠 算法 的 改造 ， 设 计 出 全 新 的 求解 0-1 背包 

问题 的 混合 蝙蝠 算法 。 有 具体 步骤 如 下 : 


At—»0, 7 »o 


0O) 局 部 搜索 。 若 rand > ri， 根 据 式 (7) 对 最 优 解 集中 选 
取 的 最 优 解 进行 局 部 搜索 , 得 出 局 部 新 解 xnew。 
d) 更 新 满意 解 。 若 rand < 4i; 且 F(xnow) <F(x,)， 更 新 当 
前 满意 解 及 其 对 应 的 目标 函数 值 。 
e) 更 新 音量 和 脉冲 发 生 率 。 通 过 式 (8) 减 小 Ai， 式 (9) 
增 大 i。 
排列 蝙 蚁 并 找到 当前 最 优 蝙蝠 x,。 
g)N_iter = N_iter + 1， 若 N_iter > N_gen， 则 结束 搜索 
输出 结果 , 否则 转 b)。 
注 : 在 初始 化 蝙蝠 种 群 、 蝙 蝠 位 置 更 新 和 反 
搜索 后 均 要 使 用 


3 ”仿真 实验 与 分 析 


为 了 验证 本 文 提 出 的 混合 蝙蝠 算法 的 有 效 性 ， 文 章 选用 两 
组 算 例 测试 其 求解 性 能 ,第 1 组 算 例 选用 文献 [23] 的 9 个 0-1 背 
包 问 题 进行 测试 ， 并 将 其 求解 结果 与 自 适应 元 胞 粒子 群 算法 
(adaptive cellular particle swarm optimization, ACPSO)P3 的 测试 
结果 进行 对 比分 析 。 第 2 组 算 例 选用 文献 [24] 的 三 个 0-1 背包 
问题 进行 测试 ， 并 与 基本 蝙蝠 算法 以 及 贪心 二 进 制 蝙蝠 算法 
(greedy binary bat algorithm，GBBA) 5 进行 比较 分 析 。 
本 文 仿真 实验 所 用 硬件 环境 为 Intel% Core™ i5-3320 CPU 
@2.60GHz，2.0 GB 内 存 ，64 位 Windows 8 操作 系统 ， 编 程 语 
言 为 MATLAB R2015a。 文献 [23] 自 适应 元 胞 粒子 群 算法 的 仿 
计算 环境 为 CorerM2 Duo CPU 为 2.93 GHz， 内 存 1.93 GB， 操 
作 系 统 Windows XP, 编程 语言 为 MAILAB 。 文献 [24] 贪心 二 


A 


算 子 局 部 
贪心 策略 进行 解 的 可 行 化 和 充分 利用 。) 


a 人 初始化 蝙蝠 种 群 。 蝙 蝠 种 群 大 小 为 *， 第 站 只 蝙蝠 的 位 置 
为 xi， 速度 为 Vis 追随 概率 prov， 反 置 概率 pnot， 音量 为 40， 脉 
冲 发 生 率 为 六， 最 大 友 代 次 数 N_gen， 并 根据 式 〈4) 初始 化 蝙 
蝠 位 置 


进 制 蝙蝠 算法 的 仿真 实验 所 用 的 硬件 环境 为 Intel @ Core TY 
Duo CPU T2400 1.83 GHz， 内 存 为 1GB， 操 作 系统 为 Windows 
XP2002， 编 程 语言 为 VC++ 6.0。 

为 体现 算法 对 比 结果 的 公平 性 ， 本 文 求解 第 1 组 算 例 的 蝙 


pb) 蝙蝠 速度 和 位 置 更 新 。 计 算得 出 最 优 蝙蝠 X, ， 运 用 式 


晤 种群 大 小 s、 最 大 迭代 次 数 N_gen、 单 个 算 例 运行 次 数 T 和文 


(5) 更 新 蝙蝠 速度 ， 利 用 遗传 算法 的 交叉 机 制 更 新 蝙蝠 位 


BE 
Xnewo 


献 [23] 的 粒子 群 算法 中 的 种 群 大 小 、 最 大 友 代 次 数 及 单个 算 例 
运行 次 数 设置 相同 , 即 s = 50，N_gen = 500, T = 30。 其 他 参 
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T 


数 遵从 算法 的 具体 特征 ， 设 置 为 : 初始 音量 4? = 0.25， 初 始 脉 
冲 发 生 率 72 = 0.5, a = y = 0.9， 追 随 概率 pror = 0.5， 反 置 概 
率 pnot = 0.2。 


万 晓 琼 ， 等 : 求 


表 1 给 出 第 1 组 9 个 测试 算 例 的 维 数 、 物 品 重量 集 w、 物 
品 价值 集 vs、 背 包 载 重量 C 以 及 已 知 最 优 值 。 


表 1 0-1 背包 问题 的 9 组 仿真 算 例 
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每 个 算 例 独立 测试 30 次 , 结果 如 表 2 所 示 。 其 中 ,实验 最 
优 值 、 实 验 平均 值 、 实 验 最 差 值 由 30 次 运行 结果 对 比 得 出 ; 总 
消耗 时 间 为 30 次 独立 实验 的 运行 总 时 间 ; 成 功 次 数 为 30 次 测 
试 中 找到 当前 最 优 解 的 次 数 ， 最 小 兴 代 次 数 、 平 均 迭 代 次 数 、 
最 大 迭代 次 数 为 找到 最 优 值 的 前 提 下 对 比 实验 结果 得 出 。 表 2 
中 ACPSO 算法 的 实验 结果 均 来 自 文献 [23]。 

由 表 2 可 看 出 , 对 于 算 例 KP9, HBA 算法 得 出 的 实验 最 优 
解 为 4987， 实 验 平均 值 为 4986.5， 优 于 ACPSO 给 出 的 已 知 最 
优 值 4986, 说 明 HBA 的 全 局 搜索 能 力 优 于 ACPSO; 对 于 算 例 
KP1~KP8, 从 成 功 次 数 可 以 看 出 ,HBA 运行 的 结果 除 算 例 KP6 
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以 外 ， 其 他 算 例 30 次 测试 均 能 达到 最 优 值 ， 而 ACPSO 只 有 
KP1、KP5、KP8 能 100% 达 到 最 优 值 ， 二 者 对 比 说 明 HBA 
优 的 鲁 棒 性 较 强 ， 且 从 迭代 次 数 的 三 个 实验 指标 可 以 看 出 ， 
HBA 能 在 较 低 的 迭代 次 数 中 达到 最 优 解 ， 除 KP5 算 例外 ， 均 
远 优 于 ACPSO, 说 明 HBA 具有 较 快 的 收敛 速度 。 综合 HBA 的 
成 功 次 数 以 及 寻 优 成 功 的 迭代 次 数 两 个 方面 分 析 , HBA 消耗 的 
总 时 间 显 然 优 于 ACPSO。 综合 所 有 的 实验 指标 可 以 看 出 , HBA 
的 寻 优 能 力 远 优 于 ACPSO，, 在 求解 0-1 背包 问题 中 具有 较 强 的 
求解 性 能 。 
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表 2 实验 结果 对 比 


千 例 。 千 法 维 数 。 当前 最 ”实验 最 。 实验 平 。 ”实验 最 ”消耗 总 ”成 功 ” 最 小 闪 ”平均 迁 ”最 大 先 
ee 优 值 。 。 优 值 均值 差 值 。 时 间 /s 次数。 代 次 数 “” 代 次 数 。 代 次 数 
ACPSO 295 285 295 0.19 30 180 6 
KP1 10 295 
HBA 295 285 295 0.63 30 1.00 1 
ACPSO 1024 1.0212e +003 1018 50.44 16 194.19 441 
KP2 20 1024 
HBA 1024 1024 1024 0.62 30 1.43 4 
ACPSO 1042 1.0393e +003 1037 56.45 14 49 235.00 474 
KP3 20 1042 
HBA 1042 1042 1042 3.29 30 27223 142 
ACPSO 4882 4. 8571e +003 4834 131.90 5 91 133.93 421 
KP4 50 4882 
HBA 4882 4882 4882 1,51 30 6.73 39 
ACPSO 15170 15170 15170 0.94 30 1.00 小 
KPS 100 15170 
HBA 15170 15170 15170 0.91 30 1.83 这 
ACPSO 26559 2.6527e +004 26525 367.52 1 45 45.00 45 
KP6 100 26559 
HBA 26559 2.6551e+004 26534 35.47 14 2 10.43 103 
ACPSO 2660 2660 2660 203.25 13 20 190.38 417 
KP7 100 2660 
HBA 2660 2660 2660 1.27 30 1 2.47 5 
ACPSO 4143 4143 4143 37,73 30 4 58.83 162 
KP8 100 4143 
HBA 4143 4143 4143 1.45 30 3 4.57 9 
ACPSO 4986 4986 4986 0.80 30 1 1.00 1 
KP9 100 4986 
HBA 4987 4.9865e+003 4986 31.17 15 这 4.73 33 
为 充分 测试 各 算法 的 求解 性 能 ， 体 现 对 比 公平 性 ， 本 文 在 。 ” 标 函 数值 的 比较 。 对比 50 次 测试 结果 得 出 最 差 值 、 平 均值 和 最 


测试 第 2 组 算 例 时 ， 混 合 蝙蝠 算法 中 初始 音量 4 = 0.005、 初 。” 优 值 ; 表 5 为 BA、GBBA 与 HBA 求解 算 例 1~3 的 收敛 性 比 
台 脉冲 发 生 率 r? = 0.75、 音 量 衰减 系数 a = 0.95、 脉 冲 发 生 率 ”” 较 ， 当 算法 在 最 大 迭代 次 数 限 制 内 达到 最 优 值 时 为 寻 优 成 功 ， 
增长 系数 y = 0.7、 单 个 算 例 测试 次 数 T = 50 均 沿 用 文献 24] 的 “记录 下 迭代 次 数 , 对 比 50 次 运行 结果 得 出 最 小 数 、 平 均 数 和 成 
设置 ， 追 随 概 率 prof、 反 置 概率 pnor 同 第 1 组 算 例 。 第 2 组 算 ”功率 。 表 4、 表 5 中 GBBA 的 实验 结果 均 来 自 文献 [24] 。 对 于 
例 选 用 文献 [24] 的 3 个 算 例 , 其 中 算 例 1、 算 例 2 和 第 1 组 算 例 算 例 3, 在 目标 函数 值 的 对 比 中 最 大 迭代 次 数 设置 为 300, 在 收 
中 的 KP1、KP2 相同 ， 有 具体 见 表 1; 算 例 3 如 表 3 所 示 。 敛 比 较 中 最 大 迭代 次 数 设置 为 500， 为 了 实验 对 比 结果 的 公平 
对 3 个 算 例 分 别 进行 测试 ， 结 果 如 表 4、 表 5 所 示 。 表 4 性， 本 文 和 文献 [24] 设置 相同 。 
为 BA、GBBA 与 HBA 求解 算 例 1~3 达到 最 大 进 代 次 数 时 的 目 
表 3 0-1 背包 问题 的 仿真 


~ 


算 例 维 数 背包 问题 己 知 最 优 值 


w = [438,754,699,587,789,912,819,347,511,287,541,784,676,198,572,914,988, 
4,355,569,144,272,531,556,741,489,321,84,194,483,205,607,399,747,118,651, 
806,9,607,121,370,999,494,743,967,718,397,589,193,369];v=[72,490,651,833,8 
83,489,359,337,267,441,70,934,467,661,220,329,440,774,595,98,424,37,807,32 
0,501,309,834,851,34,459,111,253,159,858,793,145,651,856,400,285,405,95, 
391,19,96,273,152,473,448,231];C=11258; 


算 例 3 50 16102 


表 4BA、GBBA 与 HBA 求解 算 例 1~3 的 目标 函数 值 比 较 


算 例 算法 天 杭 最 差 值 ”平均 值 ” 最 优 值 
BA 158 226. 42 295 
算 例 1 GBBA 4 300 294 294. 80 295 
HBA 294 294. 90 295 
BA 813 918. 36 995 
算 例 2 GBBA 4 300 1018 1019 1024 
HBA 1018 1023. 40 1024 
BA 7989 10126. 12 12224 
算 例 3 GBBA 15 300 16012 16053 16102 
HBA 16029 16086. 84 16102 
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表 5 BA、GBBA 与 HBA 求解 算 例 153 的 收敛 性 比较 


8 种 群 

算 例 算法 规模 
BA 

算 例 1 GBBA 15 
HBA 
BA 

算 例 2 GBBA 15 
HBA 
BA 

算 例 3 GBBA 15 
HBA 


表 4 的 最 优 值 可 以 看 出 HBA 和 GBBA 均 能 求 得 算 例 
1~3 的 最 优 解 ， 而 基本 BA 算法 仅 能 求 得 算 例 1 的 最 优 解 ， 表 
明 改 进 BA 算法 一 定 程度 上 克服 了 BA 易 陷 入 局 部 最 优 的 缺点 ， 
寻 优 能 力 均 优 于 BA。 由 表 4 的 平均 值 可 以 看 出 HBA 的 求解 结 
果 均 优 于 另外 两 种 算法 ， 且 从 表 5 的 成 功率 可 以 看 出 ， 对 于 算 
例 1 和 算 例 2，HBA 算法 均 能 100% 求 得 最 优 解 ， 求 解 算 例 3 
的 成 功率 也 优 于 其 他 两 种 算法 ， 综 合 两 个 实验 指标 可 以 看 出 
HBA 的 平均 求解 精度 和 求解 鲁 棒 性 均 优 于 BA 和 GBBA。 综合 
表 4、5 可 以 看 出 ， HBA 求解 0-1 背包 问题 的 综合 能 力 优 于 BA 
和 GBBA。 
为 了 更 直观 地 展现 出 HBA 的 求解 效果 ,图 2~4 给 出 BA 和 
HBA 求解 算 例 1~3 的 迭代 过 程 图 。 如 图 所 示 ，HBA 能 在 较 小 
的 迭代 次 数 内 求 得 最 优 值 ， 而 BA 在 进行 一 定 的 迭代 次 数 后 就 
陷入 了 局 部 最 优 ， 由 此 可 以 看 出 HBA 克服 了 BA 易 早熟 以 及 
收敛 速度 慢 的 缺点 ， 全 局 搜索 能 力 和 收敛 速度 均 远 优 于 BA。 
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图 2 求解 算 例 1 的 迭代 过 程 图 
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图 3 求解 算 例 2 的 迭代 过 程 图 
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结束 语 


本 文 针 对 0-1 背包 问题 的 具体 特征 ， 在 基本 蝙蝠 算法 的 框 


基础 上 加 入 遗传 算法 、 反 置 算 子 和 贪心 策略 ， 利 用 四 者 各 自 
点 ， 实 现 优势 互补 ， 设 计 出 一 种 全 新 的 混合 蝙蝠 算法 。 该 算 


法 克服 了 基本 蝙蝠 算法 易 陷入 局 部 最 优 和 收敛 速度 慢 的 缺点 。 
通过 对 0-1 背包 算 例 的 仿真 对 比 研究 ， 本 文 提 出 的 混合 蝙蝠 算 


和 贪心 二 进 制 蝙蝠 算法 , 是 求解 0-1 背包 问题 的 一 种 有 效 方法 。 


法 的 综合 求解 能 力 优 于 自 适 应 元 胞 粒子 群 算法 、 基 本 蝙蝠 算法 


参考 文献 : 


[1] 


[2] 


[4] 


Fayard D, Plateau G. Resolution of the 0-1 knapsack problem comparison 
of methods [J]. Mathematical Programming, 1975, 8 (1): 272-307. 

盛 红 波 ， 孙 娟 ， 孙 小 玲 . 0-1 多 项 式 背 包 问 题 的 一 种 精确 算法 [中 上海 
大 学 学 报 : 自然 科学 版 , 2006, 12 (4): 389-393. (Sheng Hongbo, Sun Juan 
Sun Xiaoling. A rigor method for solving 0-1 polynomial knapsack problem 
[J]. Joural of the Shanghai University: Natural Science, 2006, 12 (4): 389- 
393.) 

Shen Jingcheng, Shigeoka K, Ino F, et al. An out-of-core branch and bound 
method for solving the 0-1 knapsack problem on a GPU [J]. International 
Conference on Algorithms & Architectures for Parallel Processing, 2017: 
254-267. 

Elkihel M, Baz D E. An efficient dynamic programming parallel algorithm 
for the 0-1 kanpsack problem [J]. Parallel Computing-advances & Current 


Issues-the International Conference Parco, 2015: 298-305. 


录用 稿 


[5] 徐 颖 . 回溯 法 在 0-1 背包 问题 中 的 应 用 [J]. 软件 导 刊 ,2008 (12): 54-55. 
(Xu Ying. Application ofbacktracking method in 0-1 knapsack problem [J]. 
Software Guide, 2008 (12): 54-55. ) 

[6] 史 岚 , 张 义 宏 , 吕 建 辉 , 基于 绝对 贪心 和 预期 效率 的 0-1 背包 问题 优化 
四 . 计算 机 应 用 研究 , 2014, 31 (3): 684-687. (Shi Lan, Zhang Yihong, Lv 
Jianhui. Optimization algorithm of 0-1 knapsack problem based on absolute 
greedy and expected efficiency [J]. Application Research of Computers, 
2014, 31 (3): 684-687. ) 

[7] 胡 小 兵 , 黄 席 栅 . 基于 蚁 群 优化 算法 的 0-1 背包 问题 求解 [四 .系统 工 
程 学 报 , 2005, 20 (5): 520-523. (Hu Xiaobing, Huang Xiyue. Solving 0-1 
knapsack problem based on ant colony optimization algorithm [J]. Journal 
of Systems Engineering, 2005, 20 (5): 520-523. ) 

[8] 秦 玲 ,白云 ， 章 春芳 等 . 解 0-1 背包 问题 的 蚁 群 算法 [J]. 计算 机 工程 ， 
2006, 32 (6): 212-214. (Qin Ling, Bai Yun, Zhang Chunfang, et al. Ant 
colony algorithm for 0-1 knapsack problem [J]. Computer Engineering, 
2006, 32 (6): 212-214. ) 

[9] 李 若 平 ,欧阳 海滨 ,高 立 群 ， 等. 学 习 型 和 声 搜索 算法 及 其 在 0-1 背包 
问题 中 的 应 用 [J]. 控制 与 决策 , 2013, 28 (2): 205-210. (Li Ruoping, 
Ouyang Haibin, Gao Liqun, et al. Learned harmony search algorithm and its 
application to 0-1 knapsack problems [J]. Control and Decision, 2013, 28 
(2): 205-210. ) 

[10] Wang Ling, Yang Ruixin, Xu Yin, et al. An improved adaptive binary 
Harmony Search algorithm [J]. Information Sciences, 2013, 232 (5): 58-87. 

[11] Zou Dexuan, Gao Liqun, Li S, et al. Solving 0-1 knapsack problem by a 
novel global harmony search algorithm [J]. Applied Soft Computing Journal, 
2011, 11 (2): 1556-1564. 

[12] 王 莉 ， 绍 定 宏 ， 陆 金 桂 . 基于 遗传 算法 的 0//1 背包 问题 求解 [四 . 计算 
机 仿真 , 2006, 23 (3): 154-156. (Wang Li, Shao Dinghong, Lu Jingui. The 
genetic algorithm of solving 0//1 knapsack problem [J]. Computer 
Simulation, 2006, 23 (3): 154-156. ) 

[13] Caro F. A flipping local search genetic algorithm for the multidimensional 
0-1 knapsack problem [Cl]// Proc of Spanish Association Conference on 
Current Topics in Artificial Intelligence. 2005: 21-30. 

[14] Zhang Guoli, Wei Yi. An improved particle swarm optimization algorithm 


for solving 0-1 knapsack problem [J]. International Conference on Machine 


Learning & Cybernetics, 2008, 2: 915-918. 

[15] Haddar B, Khemakhem M, Rhimi H. A quantum particle swarm 
optimization for the 0—1 generalized knapsack sharing problem [J]. Natural 
Computing, 2016, 15 (1): 153-164. 

[16] 吴 虎 胜 ， 张 凤 鸣 ， 战 仁 军 ,等 . 求解 0-1 背包 问题 的 二 进 制 狼 群 算法 加 
系统 工程 与 电子 技术 , 2014, 36 (8): 1660-1667. (Wu Husheng, Zhang 
Fengming, Zhan Renjun, et al. A binary wolf pack algorithm for solving 0- 
1 knapsack problem [J]. Systems Engineering and Electronics, 2014, 36 (8): 
1660-1667. ) 

[17] Abdel-Basset M, El-Shahat D, Sangaiah A K. A modified nature inspired 
meta-heuristic whale optimization algorithm for solving 0-1 knapsack 
problem [J]. International Journal of Machine Learning & Cybernetics, 2017 
(1): 1-20. 

[18] Abdel-Basset M, Luo Qifang, Miao Fahui, et al. Solving 0-1 knapsack 
problems by binary dragonfly algorithm [Cl]// Proc of International 
Conference on Intelligent Computing. 2017: 491-502. 

[19] Cotta C, Troya J M. A hybrid genetic algorithm for the 0-1 multiple knapsack 
problem [J]. Artificial Neural Nets & Genetic Algorithms, 1998: 250-254. 

[20] Hu Zhijun, Li Rong. Ant colony optimization algorithm for the 0-1 knapsack 
problem based on genetic operators [C]// Advanced Materials Research. 
2011: 973-977. 

[21] Zhao Jiangfei, Huang Tinglei, Pang Fei, et al. Genetic algorithm based on 
greedy strategy in the 0-1 knapsack problem [Cl]// Proc of International 
Conference on Genetic & Evolutionary Computing. 2010: 105-107. 

[22] Yang Xinshe. A new metaheuristic bat-inspired algorithm [J]. Computer 
Knowledge & Technology, 2010, 284: 65-74. 

[23] 李 枝 勇 ,马良 ， 张 惠珍 . 求解 0//1 背包 问题 的 自 造 应 元 胞 粒子 群 算法 
四. 计算 机 工程 , 2014, 40 (10): 198-203. (Li Zhiyong, Ma Liang, Zhang 
Huizhen. Adaptive cellular particle Swarm algorithm for solving 0//1 
knapsack problem [J]. Computer Engineering, 2014, 40 (10): 198-203. ) 

[24] 吴 陪 聪 ， 贺 费 朝 ， 陈 迪 瑛 ， 等 . 求解 0-1 背包 问题 的 二 进 制 蝙蝠 算法 [J]. 
计算 机 工程 与 应 用 ,2015, 51 (19): 71-74. (Wu Congcong, He Yichao, Chen 
Yiying, et al. Binary bat algorithm for solving 0-1 knapsack problem [J]. 


Computer Engineering and Application, 2015, 51 (19): 71-74. ) 


